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ABSTRACT 


In  Section  I a status  report  is  given  on  the  possible 
'Verification  of  ARL's  forward  scattering  theory  with  data 
from  the  sea  surface.  The  literature  on  forward  and 
specular  scattering  has  been  thoroughly  examined.  No 
appropriate  forward  scatter  data  were  found  from  the  sea 
surface.  Due  to  the  lack  of  available  sea  surface  data, 
a forward  scatter  esqperlment  is  presently  being  proposed 
at  ARL's  Leike  Travis  Test  Station.  A brief  discussion 
and  proposed  esqperlmental  setup  are  described  herein,  along 
with  sonte  theoretical  forward  scatter  predictions.  In 
Section  II  the  acoustic  surface  duct  propagation  problem 
for  an  Irregular  boundary  is  solved  in  terms  of  a contour 
integral.  The  rou^^ess  is  Incorporated  into  the  boundary 
value  problem  by  means  of  an  effective  reflection  coeffi- 
cient. Application  is  meule  to  the  surface  duct  transmission 
problem  for  a velocity-depth  variation  given  by  the  Epstein 
profile . 
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I.  VERIFICATION  OF  THE  AEL  SCATTERIMG  MODEL 


The  search  in  the  literature  for  forward  scattering  data  to  be 
used  to  validate  ARL's  scattering  nodel  continued  and  was  completed 
this  quturter.  The  conqplete  bibliography  is  given  in  the  appendix  of 
this  report.  The  bibliography  is  not  a complete  listing  of  papers  on 
scattering  from  the  sea  surface;  theoretical  papers  cuid  those  con- 
cerned primarily  with  backscattering  or  reverberation  were  omitted. 
However,  the  list  does  represent  a closed  set  of  eurticles  and  reports 
on  forward  or  speculajp  scattering  from  rough  surfaces.  In  other  words, 
the  reports  eure  extensively  cross-referenced  within  the  bibliography 
itself,  with  few  extraneous  references  to  papers  not  included  in  the 
list. 


Very  few  of  these  reports  were  found  to  contain  pertinent 
scattering  data.  The  reports  and  papers  that  did  present  forward  or 
specular  scattering  data  were  found  to  be  unsuitable  for  comparison 
to  ARL's  theoretical  predictions  for  several  reasons.  The  most 
common  reasons  were  that  the  data  had  been  taken  from  model  surfaces 
or  that  certain  vltail  peuraffleters  (such  as  measurements  of  the  surface 
statistics)  were  not  Included  in  the  report.  The  model  surface  data 
in  some  cases  were  quite  good  but  were  rejected  because  the  ARL 
scattering  predictions  have  sJready  been  campeured,  with  good  results, 
to  the  scattering  data  taken  from  ARL's  model  surfaces. 


In  the  absence  of  appropriate  forward  scatter  data  in  the 
literature,  am  e3q>eriment  has  been  designed  amd  is  being  csirrled  out 
to  obtadn  forwaord  scatter  data  at  the  Ledce  Travis  Test  Station  of  ARL 


•> 


i 

! 


; 


ii 


Usual  forweird  scatter  experiments  involve  measuring  the  scattered 
field  over  a range  of  receiver  angles  with  a fixed  incident  grazing 
angle.  This  series  of  measurements  is  then  repeated  for  different 
incident  angles.  However,  this  type  of  experiment  is  difficult  to 
carry  out  in  an  ocean  or  lake  environment  hecause  of  the  changing 
receiver  position  and  the  large  distances  involved.  More  precise 
meas\irements  of  the  incident  and  receiver  angles  and  the  related 
distances  can  he  obtained  when  the  scattered  field  from  model  surfaces 
is  measured,  where,  for  example,  the  receiver  can  be  placed  on  a 
rotating  boom,  as  in  the  ARL  scattering  ejqjeriments.  The  difficulty  of 
measuring  these  angles  and  distances  and  of  obtaining  accurate  measure- 
ments of  the  surface  statistical  parameters  are  two  of  the  prime 
factors  limiting  the  amount  of  forward  scatter  data  that  has  been 
tedien  at  sea. 

The  Lake  Travis  experiment  has  been  designed  so  that  the  receiver 
position  remains  fixed  throughout  a single  data  nin.  The  source 
position  must  also  remain  fixed,  but  the  source  will  be  rotated  to 
change  the  incident  grazing. angle.  Thus,  it  is  only  necessary  to  loiow 
the  depth  of  the  source  and  receiver,  the  tilt  angles  of  the  source, 
and  the  source -receiver  separation  to  medce  a single  data  run.  Rotating 
the  source  will  change  the  direction  of  the  specularly  reflected 
component,  and,  hence,  the  receiver  will  assvane  several  different 
angles  with  respect  to  the  specular  direction  in  the  course  of  a single 
run.  The  geometry  of  this  experiment  is  shown  in  Dwg.  AS-71-1095* 
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In  order  to  convert  the  theoretical  formulas  in  the  ARL  model  to 
this  particular  geometry,  it  will  be  necess8Lry  to  find  the  new  relation 
between  the  Incident  grazing  angle,  p,  and  the  receiver  angle,  y. 

Also  the  changing  pathlengths  must  be  correctly  accounted  for  (the 
total  pathlength,  as  well  as  the  source  and  receiver  distances,  were 
previously  assumed  constant  for  any  combination  of  incident  euid 
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receiver  angles).  These  relations  can  be  easily  worked  out  in  terms 
of  the  distances  and  angles  given  in  Dwg.  AS -71-109^.  For  the  case 
of  specular  reflection  the  incident  grazing  angle  is  given  by 


P = tan 


-1  /ai 


(f) 


and  the  total  pathlength  is 


pathlength  (specular)  = 2 


(1) 


For  the  nonspecular  case, 

) 

P = I - * , (2) 

where  $ is  the  known  source  tilt  angle,  so 

i = ^ 

tan  p * 

and 

Finally  the  angle  between  the  specular  component  and  the  received 
component  is  given  by 

a = p - r , (4) 

8uid  the  total  pathlength  is  given  by 

pathlength  = Vh^  + (L-i  + Vh^  + ^ (5) 


4 


vhere  the  first  term  represents  the  distance  from  the  receiver  to  the 
point  of  reflection  on  the  surface,  euid  the  second  term  is  the 
distance  from  that  point  to  the  source. 

These  relations  have  been  incorporated  into  the  formulas  for 
the  scattering  coefficients  given  in  the  Final  Report  under 
Contract  N00024-69-C-1275.  The  scattering  coefficient  vas  calculated 
for  an  exponential  distribution  of  heights  on  the  svirface  and  a 
Gaussian  correlation  function,  using  Eq.  (75)  and  Eq.  (79)  of  that 
report.  The  resulting  plots  of  the  scattering  coefficient  versus 
incident  grazing  angle,  P,  are  shown  in  Dwgs.  AS-7I-IO96,  AS-71-1097# 
euad  AS-7I-IO98  for  rms  surface  heights  of  O.0875  in.,  0.284  in.  and 
0.875  in.,  respectively.  All  of  these  plotj  are  for  a single  scat- 
tering geometry  in  which  the  source  and  receiver  axe  separated  by  a 
distance  of  100  ft  and  are  both  submerged  at  a depth  of  50  ft.  Ihese 
graphs  show  that  the  theoretical  cxirves  are  reasonably  sensitive  to 
changes  in  the  rms  surface  helc^t  and  therefore  should  provide  a good 
confer Ison  with  experimental  data. 
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SCATTERING  COEFFICIENT  - dB 


TOTAL  COEFFICIENT 
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RANGE:  100  ft  DEPTH:  50  ft 
FREQUENCY:  100  kHz 
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II.  SOUND  PROPAGATION  IN  A SURFACE  DUCT  WITH  A RaiGH  BOUNDARY 

The  purpose  of  this  section  is  to  show  how  existing  contour 
integral  solutions  ceui  be  used  in  solving  the  surface  duct  propaga- 
tion problem  when  the  surface  is  irregular.  This  method  is  different 

1 2 

from  that  used  previously  by  Bucker.  * 

The  acoustic  field  is  derived  for  a cw  point  source  located 
in  an  inhomogeneous  medium  with  a rovigh  surface.  Application  is  made 
to  a class  of  surface  duct  profiles  in  which  the  velocity-depth  varia- 
tion is  almost  isovelocity  near  the  sxirface  (as  shown  in  Dwg.  AS-71-1099)* 
The  bilinear  and  Epstein  profiles  fit  this  class  of  velocity  variations 
very  well. 

A.  Green's  Function  Solution 

The  wave  equation  for  a point  source  with  angular  frequency  <o 
located  in  a layered  inhomogeneous  medium  is  given  by 

6(r-r  ) e"^“*  , (6) 

C^?t  ° ° 

where  P represents  the  pressure,  C(z)  is  the  sound  velocity  (variable 
in  the  depth  coordinate  Z),  and  6(r-f^)  is  the  three-dimensional 
Dirac  delta  function.  The  strength  of  the  source  is  designated  by 
If  cylindrical  coordinates  (r,B,Z)  are  assianed  to  have  azimuthal 
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s^BMtry  and  a tine  factor  of  exp(-iait)  is  suppressed,  the  following 
equation  Is  obtained: 


id/  dP\ 


d^p 


k^(Z)P 


-2Q^6(r)B(Z-Z^) 


(T) 


where  P now  represents  the  time  Independent  pressure.  The  wave 

nvmber  k(z)  is  defined  as  a/C(Z)  and  the  source  is  located  at 

Z » Z and  r = 0,  as  indicated  in  Dwg.  AS-71-1099.  The  vertical  depth 
o 

coordinate  Z varies  from  Z^^  * Z S <»,  and  the  range  coordinate  r varies 
from  0 < r < ».  Boundary  conditions  inposed  on  Eq.  (7)  sure  that 


(a)  P must  satisfy  a radiation  condition  for  r 

(h)  effective  upgolng  and  downgoing  reflection  coefficients 

sure  specified  at  the  surface  Z = Z.> 
dp  ^ 

(c)  P and  ^ he  continuous  across  any  discontinuities  in  the 
velocity-depth  profile . 

The  statistical  nature  of  the  surface  is  Introduced  through  one  of 
the  effective  reflection  coefficients. 


In  the  Green's  function  approach,  Eq.  (7)  is  separated  into  the 
following  forms: 


[5  (’■  h)  * v] 

^ + k^(Z)-\^l  Gg(Z,Z^,-\^)  = -<l^6(Z-Z^)  , (9) 
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where  is  the  separation  constant.  By  the  resolvent  Green's 
function  technique,  the  solution  of  Eq.  (7)  is  given  by  a complex 
convolution  of  Gj^  and  G^, 

P(r,Z,Z„)  . ^ / Sl(r,Xi)  , (10) 

where  the  contour  C separates  the  singularities  of  the  Green's  functions. 

The  solution  of  Eq.  (8)  which  satisfies  the  required  boundary 
conditions  is  given  by 

Gi(r,\^)  = ijtH^(lr)  0 < \ , (ll) 

where  H^(k  rs)  is  the  Hankel  function  of  the  first  kind,  and  where 


I = k sin  6 = k s , 
o o' 


0 is  the  angle  of  incidence,  and  k^  is  the  wave  number  at  the  surface 
Z = Z . The  solution  of  Eq.  (9)  that  satisfies  the  required  boundary 
conditions  is  given  by 


G(Z,Z^,s)  - 


.1 

2k  (l-s  ) 
o 


G(Z,Z  ,s)  = 

o o 

2k^(l-s^) 


. fp(Z)f-(Z) 

-TT75  1 - VTT 
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In  order  to  define  the  functions  in  Eq.  (13)>  Dwg.  AS-71-1099  is 
separated  as  shown  in  the  following  sketch. 


The  idea  is  to  divide  an  inhomogeneous  layer  into  two  profiles  where 
each  has  a homogeneous  half-space.  In  the  present  application, 
region  I of  sketch  (c)  is  also  homogeneous.  is  defined  as 

the  reflection  coefficient  of  a plane  wave  incident  from  the  homo- 
geneous medium  onto  the  lower  Inhomogeneous  half- space  [sketch  (b)]. 
The  function  fj^(Z)  is  the  field  in  the  inhomogeneous  half-space 
II -b,  when  the  incident  wave  has  unit  amplitude.  Vg(0)  and  fg(Z) 
have  an  analogous  meaning  in  sketch  (c).  At  the  source  = 1+Vj^ 

and  fg(Z^)  = l+Vg. 

When  Eq.  (ll)  and  Eq.  (13)  are  used  in  Eq.  (lO),  the  solution  for 
Eq.  (7)  is  obtained, 


-ik  a / f (Z^)f  (Z)  , ^ 

P(r,Z,Z„)  . V“>  ^ h ^ ^ 

P(r,Z,Z„)  . Lv'V'  V“)  ^ ^ < 


- 


where  s = sin  6 and  c = cos  6. 
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For  the  case  where  is  a lossy  boxindary,  then,  from  classical 


scattering  theory  the  form  for  Vg  can  he  given  as 


W^Q)  = exp^-2(k^a  cos 


where  a is  the  rms  height  value  for  a rough  ocean  surface.  If  a 
peaJk  to  trough  height  h is  known,  then  a good  approximation  for  a is 
given  hy  the  equation 


The  wave  height  h can  be  expressed  as  a fvinction  of  the  wind  speed  W 
(in  knots)  by  the  following  relation  (derived  by  Schulkin); 

h = 0.00276  , (17) 

where  h is  given  in  feet. 

Equation  (l4)  will  now  be  applied  to  the  case  of  a rough  surface 
bounding  an  Epstein  profile.  This  profile  is  used  instead  of  the 
bilinear  because  it  does  not  have  suiy  velocity-depth  discontinuities. 

B.  Epstein  Profile  with  a Rough  Surface 


The  Epstein  profile  is  a four-parameter  velocity  function  which 
has  had  many  applications.  Since  Eckart  as  well  as  Rosen  suid  Morse 
examined  this  function  about  the  same  time  that  Epstein  did, it  is  some- 
times known  by  other  names.  Specifically,  the  velocity-depth  profile 
known  as  the  Epstein  function  is  given  by 


C (Z) 


= A sech^  5 ^ ® tanh  5 ^ 


vhere  the  qiiantities  A,  B,  D,  and  H are  profile  parameters.  Bucher 
and  Morris^  have  given  a normal  mode  treatment  for  an  Epstein  surface 
duct  with  a pressure  release  boundary.  The  profile  shown  in 
Dwg.  AS-7I-IIOO  has  been  taken,  with  permission,  from  the  paper  by 
Bucker  and  Morris. 


Epstein  chose  this  function  because  it  allows  Eq.  (9)  to  be 
solved  in  terms  of  hypergeometric  functions.  Specifically,  the 
traveling -waves  solutions  to  the  homogeneous  Eq.  (9)  are  given  by 

n^(Z)  . fJi-v  . i|!i.  V * 1 ♦ M *)  (19) 

and 

nJiZ)  - (1)^^  Fg^l-v  + ^,  V + 1 + p;  1 - xj  , 

(20) 


where  F(x)  li  the  hypergeometric  function,  and 


K = , 

= h(^-  kf.) 

, 

V . i . , 

C = (D  + B)"^/^  = C(Z  = -ho) 
00 

C = Z/H  , 

C ^ = (D-  B)“^/^  = C(Z  = -« 

X = (l-f  tanh  C)/2  , 

C = (D  + A)“^/^  = C(Z  = 0) 
0 

(21) 
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WATER  DEPTH  - yd 


The  following  limiting  forms  are 


Urn  n (Z)  ~ (H/2)^/^  exp(\?) 
z -♦-» 

lim  n„(z)  ~ (H/2)^/^  exp(-MC) 


z -»+» 

2 2 2 2 
Also  when  k_  > | and  k > | , then 


/ 2 2 / 2 

=4^-0 

and  (2! 

/ 2 2\^/^  / 2 2\^/^ 

In  order  to  find  the  pressvire  field  in  an  Epstein  surface  duct, 
the  appropriate  Green’s  function  must  be  specified.  Equation  (13) 
will  now  be  given  in  a more  general  form. 


The  following  diagram  is  necessary  for  the  definition  of  the 
Green's  function. 


When  ^ Zg,  the  Green's  function  is  given  by 

f n-  (Z,  )n  JZ  J If  n (Z  JngCZg)’ 

K<^<)  * ni^)  -|h<V  * 

■ r-^ — ..■(iijn^u-s-ri ^ 

W(ng,n^)  1 - R^(Z^)Rg(Z2^  S^Tzpn^Tz^ 
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where  emd  denote  the  smaller  or  larger,  respectively,  of  the 
variables  Z and  Z^.  V{n^,n^)  represents  the  Wronsklan  of  ®nd 
nj^(Z)  must  be  an  outgoing  wave  for  Z and  n^(Z)  must  be  an 

outgoing  wave  for  Z The  reflection  coefficients  are  defined  as 


. 


R,(Z,)  = . 


R,(Z.)  . 


(25) 


where,  by  analytic  continuation,  the  functions  n^  and  n^  are  written  as 

+ 

(26) 


nj^  = n^^  + n^ 


ng  = ng  + ng 


(That  is,  n^  «und  1^2  ere  decomposed  into  upgolng  and  downgoing  waves.) 


If  the  region  between  Z^  «md  Z2  is  taken  as  honogeneous  and  thin, 
then  Eq.  (24)  can  be  slnpllfled.  In  the  case  of  the  Epstein  profile, 
Z^  and  Z2  ere  taken  in  the  region  where  the  velocity  variation  is 
negligible . 


When  the  reflection  coefficient  R2  is  referred  to  the  level 
Z = Z^f  then  the  eigenvalue  equation  can  be  written  as 

P^(Zj^)R2(Z^)  = 1 . (27) 

and  R2  are  generally  specified  in  terms  of  impedance  conditions 
at  the  boundaries.  However,  in  the  present  treatment  R^  and  Rg  will 
not  be  functions  of  the  boundary  impedances. 
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Equation  (24)  will  nov  be  specialized  to  the  Epstein  profile 
with  a rough  surface.  The  function  n^  and  for  the  Epstein  profile 
are  those  already  given  by  Eqs.  (19)  and  (20).  For  the  reflection 
coefficient  Eq.  (15)  will  be  used,  thus 


Rj^(Zi)  = ej5»^-2(k^  a cos  8)^)  , 


where 


k 3 k 
o 


In  order  to  obtain  Rg(Zj^),  the  analytic  continuation  for  ng(Z)  is  used 
to  give 


r(-\)r(v  + ^)r(i  -v  + ' ' 

R2(Zi)  f 

^ ^ r(\)r  / V + ^ ^ 

which  is  valid  for  2|  ^ | » 1.  The  Wronskian  of  Ug  and  n^^  has  already 
been  given  in  reference  4. 


When  the  period  equation  1 - R^Rg  = 0 is  solved  for  the  complex 
eigenvalues  {,  then  the  residues  yield  the  normal  modes.  Since  the 
eigenvalues  are  complex,  the  modes  are  leaky  (individually  they  do  not 
satisfy  the  radiation  condition).  The  following  method  is  being  used 
to  calculate  the  complex  eigenvalues. 

First  the  period  equation  is  written  as 


R-R«  exp/-2ik  Z,  cos  8 ) = e 
i.  c \ o i.  m/ 


-2ijnn 


19 


where 

\ = l\l  exp(i<jJ  , 

Rg  = 1^2!  ®*P(^*2)  • 

The  complex  eigenvalues  ^ = Rq  presently  helng 

conqputed  by  means  of  a complex  gamma  function  routine  which  confutes 
valvjes  for  Re  in  T(Z)  and  Im  in  P(z).  This  Is  fortunate,  since 
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Re  fn  r(z)  = in  |r(z)| 


Im  in  r(Z)  = arg  r(Z) 


(Z  in  the  ganna  functions  represents  a complex  number).  These 
complex  roots  ■will  then  be  used  in  calciilating  the  residues.  It  is 
expected  that  this  method  will  provide  a way  to  incorporate  the  effects 
of  a rou{^  boundary  on  surface  duct  propagation. 
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In  Section  I a status  report  is  given  on  the  possible  verifica- 
tion of  ARL's  forward  scattering  theory  with  data  from  the  sea 
surface.  The  literature  on  forward  and  specular  scattering  has 
been  thoroughly  examined.  No  appropriate  forward  scatter  data 
were  fovmd  from  the  sea  surface.  Due  to  the  lack  of  available 
sea  surface  data,  a forward  scatter  experiment  is  presently  being 
proposed  at  ARL's  Lake  Travis  Test  Station.  A brief  discussion 
and  proposed  experimental  setup  are  described  herein,  along  with 
some  theoretical  forward  scatter  predictions.  In  Section  II  the 
acoustic  surface  duct  propagation  problem  for  an  irreguleu*  bound- 
ary is  solved  in  terms  of  a contour  integral.  The  roughness  is 
incorporated  into  the  boundary  value  problem  by  means  of  an  effec- 
tive reflection  coefficient.  Application  is  made  to  the  surface 
duct  transmission  problem  for  a velocity-depth  variation  given  by 
the  Epstein  profile.  (U-FCUO) 


